Non-equilibrium transport at a dissipative quantum phase transition 
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We investigate the non-equilibrium transport near a quantum phase transition in a generic and 
relatively simple model, the dissipative resonant level model, that has many applications for nanosys- 
tems. We formulate a rigorous mapping and apply a controlled frequency-dependent renormalization 
group approach to compute the non-equilibrium current in the presence of a finite bias voltage V 
and a finite temperature T. For y ^ 0, we find that the conductance has its well-known equilibrium 
form, while it displays a distinct non-equilibrium profile at finite voltage. 
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In recent years, quantum phase transitions (QPTs) 
1, % have attracted much attention at the nanoscale 
3, 5 H, H, 0, Hi- A finite bias voltage appUcd across 
a nanosystem is expected to smear out the equilibrium 
transition, but the current-induced decoherence might 
act quite differently as compared to thermal decoher- 
ence at finite temperature T, resulting in exotic behavior 
near the transition. Non-equilibrium effects at a quan- 
tum phase transition appear as an emerging field both in 
experimental and theoretical condensed matter physics 
S l3j Hi 121 ■ Quantum impurity systems out of equi- 
librium are also extensively studied theoretically 13;]. In 
this Letter, we aim to answer several fundamental ques- 
tions related to the non-equilibrium transport in quan- 
tum dot settings, such as what is the scaling behavior 
of the conductance at zero temperature and finite bias 
voltage near the transition. 

For this purpose, we employ a typical nano-model com- 
prising a dissipative resonant level (quantum dot). In 
this model, the QPT separating the conducting and in- 
sulating phase for the level is solely driven by dissipation, 
which can be modeled by a bosonic bath. Dissipation- 
driven QPTs have been addressed theoretically and ex- 
perimentally in various systems, such as: Josephson junc- 
tion arrays 14] , superconducting thin films [15], super- 
conducting qubits [l6| and biological systems 17]. 

Our Hamiltonian takes the precise form: 



H= ^ {e{k) - ^li)cl^ck^ + Ucl^d + h.c. (1) 

fe,i=l,2 



where ti is the (real) hopping amplitude between the lead 
i and the quantum dot, Cki and d are electron operators 
for the (Fermi-liquid type) leads and the quantum dot, 
respectively, /i^ = ±y/2 is the chemical potential applied 
on the lead i, while the dot level is at zero chemical po- 
tential. To simplify the discussion, we assume that the 
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G(ac,T) ~ 1/ In^ (T/D„) 
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FIG. 1: Schematic phase diagram as a function of tempera- 
ture T, dissipation strength a, and bias voltage V . 

electron spins have been polarized through the applica- 
tion of a strong magnetic field. Here, 6q are the boson 
operators of the dissipative bath, that is governed by an 
ohmic spectral density [^: J{u)) = X^5{uj — ujr) = acu. 
We use units in which h = 1 and electric charge e = 1. 

In equilibrium (V — 0), such a dissipative system 
comprising several leads maps onto the anisotropic one- 
channel Kondo model; we denote by t the hopping am- 
plitude between the (effective) lead and the level. We 
introduce the dimensionless transverse Kondo coupling 
which is proportional to t (the exact prefactor is 
given in Refs. [J d, [B]) and the longitudinal coupling 
gi'^^ oc 1 — ^/a 0, H, @1 • The model exhibits a Kosterlitz- 
Thouless (KT) QPT from a delocalized (Kondo screened) 

(e) 

gz > 0, with a large conductance, 
l/h (e ~ 1 and h = 2'Kh = 2tt), to a localized 
(local moment) phase for g^'' -I- gi'^^ < 0, with a small 
conductance, as the dissipation strength is increased (see 
Fig. 1). For —> 0, the KT transition occurs at 
c^c = 1. As a —!■ etc, the Kondo temperature Tk obeys 
In TV cx l/(a - ac) [3]. 
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In equilibrium, the scaling functions g'^f^{T) and 
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gz^\T) at the quantum critical point are obtained via the 
renormaUzation-group (RG) equations of the anisotropic 
Kondo model: ffft(^) = -Mr{T) = (2 In (T^^/T))-!; 
hereafter, we introduce the energy scale To = Dqc^^ i^9±) ^ 
with Dq being the ultraviolet cutoff and we set the Boltz- 
mann constant ks — 1- Having in mind a quantum dot 
at resonance, Dq = min((5e, ojc), with 6e being the level 
spacing on the dot and ujc the cut-off of the bosonic bath; 
-Do is of the order of a few Kelvins. At the KT quantum 
phase transition, the conductance drops abruptly 0]: 



Geq{ac,T <t: Do) cx 



(2) 



Below, we analyze the nonequilibrium {V ^ 0) transport 
at the phase transition and in the localized phase. 

First, we envision a non-equilibrium mapping revealing 
that the leads are controlled by distinct chemical poten- 
tials. Through similar bosonization and refermionization 
procedures as in equilibrium [sl, 0j S 0, our model is 
mapped onto an anisotropic Kondo model with the effec- 
tive (Fermi-liquid) left (L) and right lead {R ) (isj : 

Hk = t*^*^ ~ Ai7]47<TCfe7<T (3) 

k,-i=L,R,a=1,l 

+ {Ji-s+RS- + J1-4lS- + h.c.) + ^ J,4^5^ 

where c[,^^^^^ is the electron operator of the effective 
lead L{R), with a the spin quantum number, S'^ = d\ 
— d, and S'^ = Q — f/2 where Q — d^d describes 
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the charge occupancy of the level. Additionally, s, 
J2a.s.k.k' ^/'^^ija'^asCk'iSS are the spin-flip operators be- 
tween the effective leads 7 and /3, cx ix(2) embody 

the transverse Kondo couplings, Jz oc 1/2(1 — l/\/2a*), 
and = ±-j^l/(2a*), where l/a* = 1 + a. It should 
be noted that this mapping is exact near the phase tran- 
sition where a ^ 1 or a* ^ 1/2, and thus /i^ = ztV/2. 

From the mapping, N1 — N2 = {Nl — Nn) , where Ni = 
Sfei c^i^fei represents the charge in lead i — 1,2, whereas 
= cl^^Ck-ycr represents the charge in the effective 
lead ^ = L,R. This allows us to check that the averaged 
currents within the Keldysh formalism are the same in 
the original and in the effective Kondo model. Thus, the 
current / can be computed from the Kondo model. 

The poor-man scaling equations of Anderson are gener- 
alized to nonequilibrium RG equations by including the 
frequency dependence of the Kondo couplings and the 
decoherence due to the steady-state current at finite bias 
voltage • For the sake of clarity, we assume that the 
resonant level (quantum dot) is symmetrically coupled 
to the right and to the left lead, ti — t2- The dimen- 
sionless Kondo couplings then have the extra symme- 
try {u is the frequency): = 9±{z){-(^) where 



9±{z) = -^(0) J_L(2) with iV(0) being the density of states 
per spin of the conduction electrons. We obtain [foj : 
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where = 0(Z) — \lo ^ zFj), D < Dq is the running 
cutoff, and F is the decoherence (dephasing) rate at finite 
bias which cuts off the RG flow [19| . The configurations 
of the system out of equilibrium are not true eigenstates, 
but acquire a finite lifetime. The spectral function of the 
fermion on the level is peaked aXuj — ±y/2, and therefore 
we have 5_l(z)(w) ~ 5_l(z)(±F/2) on the right hand side 
of Eq. Q . Other Kondo couphngs are not generated. 

In the Kondo model, F corresponds to the relaxation 
rate due to spin flip processes (which are charge flips in 
the original model). From Ref. |19|], we identify: 

r = ^Xl '^'^9l{^)U-ti^ii - fuj-t^y) (5) 
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where f^j is the Fermi function. Here, 7 = 7' for the g^i^) 
terms while 7 7^ 7' for the g']_{(^) terms with 7, 7' being 
L or R. We have introduced the occupation numbers 
for up and down spins satisfying nij + = 1 and 
Sz = l/2(n| — ni). In the delocalized phase, we get = 
nj^ — 1/2 in agreement with the quantum Boltzmann 
equation [l9j : at the phase transition we can use that 
g±i^) = ~gz{^) and that X^o- ""o^ ~ ^'^'^ finally in 
the localized phase < — g^, n„ satisfies \Sz\ — > 1/2 

Following the scheme of Ref. [19|, we solve Eqs. (jl]) 
and (O self-consistently. First, we compute g±(z){^ — 
ztV/2) for a given cutoff D. Second, we substitute the 
solutions back into the RG equations to get the general 
solutions for gj_(^z)i'-^) at finite V, and then extract the 
solutions in the limit D ^ 0. When the cutoff D is 
lowered, the RG flows are not cutoff by V but continue 
to flow for F < D < y until they are stopped for D <r. 

At the KT transition, we both numerically and ana- 
lytically solve Eqs. (g]) and ^ (in the limit of £> ^ 0): 



ff_L,cr(w) 



1 



= Y,^i\^-^V/2\-V) 

P 4 In L|^_/3y/2| 

+ Q(y- \uj-(3V/2\) X 
1 



(6) 



ln[r|,/T/max(|c.-/3l-/2|,F)] 41n^ 



The solutions at the transition (denoted g±^cr and gz,cr) 
are shown in Fig. [51 Since g^.cr{^) decreases under the 
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FIG. 2: g±^cr{i^) = —52,^(1^) at a = Oc; the bare couplings 
are g± = -5^ = 0.1. We have set V = 0.066L>o, 0.43Do, 
0.7573o and 0.99Do where Do = 1 for aU the figures. The 
arrows give the values of g±^cr{i^ ~ 0) at these bias voltages. 




G(0'c,V) by numerics 
G(Cc,V) by equilibrium scaling 
— G(ac,V)byEq. (6) 
G(ac,V)byEq. (8) 



FIG. 3: Nonequilibrium conductance at the KT transition. 
Go is the equilibrium conductance at the transition for T = 
Dq: Go — Geq[a.c,T = Dq) = O.OOSvr with the bare couplings 
g± = —gz = 0.1. Again, we set the charge e = h = 1. 



RG scheme, the effect of the decoherence leads to min- 
ima; the couplings are severely suppressed at the points 
uj = i-j. We also check that g±^cr{i^) = —gz,cr{i^)- 

From the Keldysh calculation up to second order in the 
tunneling amplitudes, the current reads: 
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duj^'^4:g±{ujfna x 

(7 

/c^-Mi(l - /c^-mr) - (L ^ R). 



(7) 



At T = 0, it simpHfies as / = f J^y^2'^uigj_{uj). Then, 
we numerically evaluate the nonequilibrium current. The 
conductance is obtained from G{V) = dl/dV. The T — 
results at the KT transition are shown in Fig. [3] First, 
it is instructive to compare the non-equilibrium current 
at the transition to the approximate expression: 



^(J[5^,„,(c. = 0)f) (8) 
^{{l~^)[g^M^^V/2f), 



where g±^cri^ = 0) 
g^^^riuj = V/2) w 1 



, and 



ln(2T^/y2) - 4ln(TD/V)^ 

We have treated f;_L. cr(w)^ 
within the interval —V/2 < oj < V/2 as a semi-ellipse. 

As demonstrated in Fig. [3l the conductance G{V) ob- 
tained via the approximation in Eq. ([8|) fits very well 
with that obtained numerically over a whole range of 
< V < Dq. In the low-bias V —f (equihbrium) 

limit, since g±^cr{^ 0) « 5j^cr(^ — V) 't^ 1, we have 

I{ac,V) « ^ [9±IAT = ^))^ therefore the scaling of 
G(ac, V) is reminiscent of the equilibrium expression in 

Eq. ©, G{a.,V) « f [g^ll^iT = V))' = f i„.(^^/v-) - 
This agreement between equilibrium and nonequilibrium 
conductance at low V persists up to a crossover scale 
V ~ O.OlDo (determined for the parameters used in Fig. 



3). At larger biases, the conductance shows a unique 
nonequilibrium profile; see Eq. ([5]). We find an excellent 
agreement of the nonequilibrium conductance obtained 
by three different ways — pure numerics, analytical so- 
lution Eq. ([5]) and the approximation in Eq. ([S]). 

The distinct nonequilibrium scaling behavior seen here 
is in fact closely tied to the non-trivial (non-linear) V 
dependence of the decoherence rate r{V). In particular, 
at the KT transition, we find that F ^ i/ is a highly 
non-linear function in V, resulting in the deviation of 
the nonequilibrium scaling from that in equilibrium. For 
large bias voltages V — s- Dq^ since g±^cr{'-^) approaches 
its bare value g±^ the nonequilibrium conductance in- 
creases rapidly and reaches G{ac, V) ^ Gq = ^g\- The 
nonequilibrium conductance is smaller than the equilib- 
rium one, G{ac,V) < Geq{ac,T = V), since g±(Lu — 
±.V/2) < g±{uj — 0). Additionally, in the delocalized 
phase for V » Tk > 0, the RG flow of g± is suppressed 
by the decoherence rate, and G cx l/ln^(y/TR-) [l9|. 

In the localized phase, the equilibrium RG equations 
of the effective Kondo model can be solved analytically, 



resulting in G\^1{T) = f 



where 



2cg^{c+\g,\) 



(C+I.9.IF 



9l& 



T 



2c 



(9) 



with c = ^gf 



9l- 



We introduce the energy scale 



T* = Doe^'^/v^Psi (which vanishes at the KT tran- 
sition) such that g'fl^ciT) cx (T/T*f'' for T 0, lead- 
ing to g|;;^(T) (X (T/T*)^". For very small bias voltages 
F ^ 0, we find that the conductance reduces to the equi- 
librium scahng: G(y) G'fJ^{T = V) (x {V/T*Y'' (see 
Fig. H (a) and Fig. 4 (b)). For g^^ioc < \gz,ioc\ and 
a* 1/2, we get that the exponent 4c « 2a* — 1, in 
perfect agreement with that obtained in equilibrium at 
low temperatures: G{T) cx T^"*"! 0. At higher bias 
voltages O.OIDq <V < Dq, the conductance now follows 
a unique nonequilibrium form (consult Fig. [D^c)). 
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FIG. 4; Conductance in the localized phase (in units of 
Ti/h). (a) G{V) at low bias follows the equilibrium scaling 
(dashed lines), (b) The conductance G{V)/Gc is a func- 
tion of V/T* where we have defined Gc = G{ctc, V) and 
T* = Doe"''/V9?-9i. (c) At large bias vohages V, the 
nonequilibrium conductance GiV) (solid lines) is distinct 
from the equilibrium form (dashed lines). The dot-dashed 
lines stem from an analytical approximation via Eq. Q. 
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FIG. 5: Scaling of the conductance at the KT transition (same 
unit as in Fig. (a). For V > T, the conductance follows 
the nonequilibrium scaling G{ac,V). (b). For V <T, now 
the conductance follows the equilibrium scaling G{oic,T). 



We have also analyzed the finite temperature profile of 
the nonequilibrium conductance at the transition. We 
distinguish two different behaviors. For V > T, the 
conductance G{V, T) follows the nonequilibrium form at 
T (see Fig. [^a)), while for F < T it follows the 
{V = 0) finite-temperature expression (see Fig. (H^b)). 
These two scaling behaviors have a crossover aX V — T . 

In summary, we have investigated the nonequilibrium 
transport at a QPT using a standard nano-model, the 
dissipative resonant level. We have used an exact map- 
ping and applied a controlled frequency-dependent renor- 
malization group approach to compute the current. For 
y ^ 0, the conductance G follows the equilibrium be- 
havior; by increasing V ^ the frequency-dependence of the 
couplings begins to play an important role and therefore 
we systematically find very distinct scalings. We have 
also analyzed the finite temperature profile of G{V, T) 
at the transition and identified two distinct behaviors at 
V > T and V < T. Finally, our results have a direct 
experimental relevance for dissipative two-level systems. 

We are grateful to D. Goldhaber-Gordon and G. 
Zarand for stimulating discussions, and to R.T. Chang 
and K.V.P. Lata for technical support. We also acknowl- 
edge the generous support from the NSC grant No.95- 
2112-M-009-049-MY3, the MOE-ATU program, the 
NCTS of Taiwan, R.O.C. (C.H.C.), the Department of 
Energy in USA under the contract DE-FG02-08ER46541 
(K.L.H.), the DFG via SFB 608 and SFB/TR-12 (M.V.), 
and the DFG-Center for Functional Nanostructures, 
C.F.N. (P.W.). C.H.C. has also benefitted from the vis- 



[1] S. Sachdev, Quantum Phase Transitions, Cambridge Uni- 
versity Press (2000). 

[2] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, 
Rev. Mod. Phys. 69, 315 (1987). 

[3] K. Le Hur, Phys. Rev. Lett. 92, 196804 (2004); M.-R. 
Li, K. Le Hur, and W. Hofstetter, Phys. Rev. Lett. 95, 
086406 (2005). 

[4] K. Le Hur and M.-R. Li, Phys. Rev. B 72, 073305 (2005). 
[5] P. Cedraschi and M. Biittiker, Annals of Physics (NY) 
289, 1 (2001). 

[6] A. Furusaki and K. A. Matveev, Phys. Rev. Lett. 88, 

226404 (2002). 
[7] L. Borda et al. , |cond -mat/0602019 
[8] G. Zarand et al, Phys. Rev. Lett. 97, 166802 (2006). 
[9] D. E. Feldman, Phys. Rev. Lett., 95, 177201 (2005). 
[10] A. Mitra et al, Phys. Rev. Lett., 97, 236808 (2006). 
[11] R. M. Potok et al, Na ture 447, 167 (20 07). 
[12] S. Kirchner and Q. Si, larXiv:0805.3717l 
[13] See, e.g., P. Mehta and N. Andrei, Phys. Rev. Lett. 96, 
216802 (2006); S. Kehrein, Phys. Rev. Lett. 95, 056602 
(2005); E. Boulat, H. Saleur, and P. Schmitteckert, Phys. 
Rev. Lett. 101, 140601 (2008); B. Doyon, Phys. Rev. 
Lett. 99, 076806 (2007). 
[14] G. Refael, E. Demler, Y. Oreg, and D. S. Fisher, Phys. 

Rev. B 75, 014522 (2007). 
[15] N. Mason and A. Kapitulnik, Phy. Rev. B 65, 220505 
(2002). 

[16] A. J. Rimberg et al, Phys. Rev. Lett. 78, 2632 (1997). 

[17] J. Gilmore and R. McKenzie, J. Phys. C. 11, 2965 (1999). 

[18] In the delocalized phase, it would be more judicious to 
redefine ^ — > {t\C\ + t^c^)! \/t\ -I- t\, and then bosonize 
and refermionize to get the fermions with spins 

[19] A. Rosch et al., Phys. Rev. Lett. 90, 076804 (2003); J. 
Phys. Soc. Jpn. 74, 118 (2005). 



